Origin of Ferroelastic Domains in Free-Standing Single Crystal Ferroelectric Films 
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The origin of the unusual 90°ferroelectric / ferroelastic domains, consistently observed in recent 
studies on meso and nanoscale free-standing single crystals of BaTiC^ [Schilling et al., Physical 
Review B, 74, 024115 (2006); Schilling et al., Nano Letters, 7, 3787 (2007)], has been considered. 
A model has been developed which postulates that the domains form as a response to elastic stress 
induced by a surface layer which does not undergo the paraelectric-ferroelectric, cubic-tetragonal 
phase transition. This model was found to accurately account for the changes in domain periodicity 
as a function of size that had been observed experimentally. The physical origin of the surface 
layer might readily be associated with patterning damage, seen in experiment; however, when all 
evidence of physical damage is removed from the BaTiOs surfaces by thermal annealing, the domain 
configuration remains practically unchanged. This suggests a more intrinsic origin, such as the 
increased importance of surface tension at small dimensions. The effect of surface tension is also 
shown to be proportional to the difference in hardness between the surface and the interior of the 
ferroelectric. The present model for surface tension induced twinning should also be relevant for 
finely grained or core-shell structured ceramics. 

PACS numbers: 77.80.Bh, 77.55.+f, 77.80.Dj 



Introduction 

In 1935 Landau and Lifshitz predicted the appearance 
of periodic, thermodynamically stable, domains, with op- 
positely oriented magnetic moments, in ferromagnetic 
crystals. The existence of these domains minimized the 
energy of the depolarizing field, caused by the abrupt 
discontinuity of the spontaneous magnetization at the 
sample surface [l|. Additional consideration of the as- 
sociated domain wall energies allowed them to predict 
a square root dependence of the domain period on the 
sample thickness. This square root relationship is often 
referred to as the Kittel law because of its independent 
formulation in 1946 by Kittel [2|]. The law is also valid 
for 180°domains in ferroelectric crystals 0,0] , where the 
unfavourable depolarizing electric field is provided by the 
abrupt polar discontinuity at the surface. The Landau- 
Lifshitz-Kittel theory was later extended by Roytburd Q 
to describe the behaviour of ferroelastic domains formed 
as a result of substrate clamping effects in thin film sys- 
tems. 

In previous studies we have shown that the Kittel law 
works perfectly for a wide class of ferroic materials (fer- 
romagnetic, ferroelectric or ferroelastic) over six orders 
of magnitude in film thickness @, @] , and that it can be 
intuitively expressed in terms of the domain wall thick- 
ness [8, 9, 10]. Moreover, we have demonstrated that the 
Kittel approach can be extended to three dimensional 
structures EH, EH, to ferroelectric superlattices @ 
and to multiferroic materials 14 , THj] . Much of the exper- 



imental work has been in association with observations of 



periodic 90°ferroelectric-elastic domains that have been 
consistently observed in free-standing single crystal thin 
films (see Fig. la) and nanowires of BaTiOs. In all our ex- 
periments, the size of the 90°domains as a function of size 
is indeed found to be well described by a Kittel- Roytburd 
formalism. 

At first glance, however, the very existence of such do- 
mains in our free-standing samples is quite surprising. 
Ferroelastic domains normally appear in response to an 
external stress (such as that imposed by clamping to a 
rigid substrate, for example) which forces the sample to- 
wards shape preservation in the clamped directions. The 
domain configuration is such that the macroscopic shape 
difference between the paraphase and the ferrophase is 
minimized, while the domain size responds to an equi- 
librium between domain energy and domain wall energy 
However, our BaTiOs lamella are free standing sin- 
gle crystals and therefore free from epitaxially-induced 
stress, so they have no Roytburd-like interface-induced 
elastic driving force. The question, then, is the following: 
if there is no external stress, what causes the appearance 
of the ferroelastic domains? 

In the present work we explain the appearance of 
the self-organized domain patterns in free-standing nano- 
samples of BaTiOs by assuming that the driving stress 
is provided by an encapsulating surface layer. In much 
of our experimental work, this encapsulation layer could 
easily be associated with the surface damage caused by 
focused ion beam milling. However, we have also ob- 
served here that 90°domains persist, with only slightly 
altered periodicities, even when surface damage has been 
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FIG. 1: (Color Online) (a) Scanning Transmission Elec- 
tron Microscopy Image of periodic domain structure in 
free-standing single crystal lamela of BaTiOs. (b) Sketch 
of our model system, comprising a ferroelectric layer 
(FL) sandwitched between two surface encapsulation lay- 
ers (SL) which do not undergo the ferroelastic transforma- 
tion. The stress imposed by these untransformed "dead" 
layers onto the ferroelectric layer induces the appearance of 
90°ferroelectric / ferroelastic domains 



repaired by thermal annealing; this suggests that the 
strain effects may in fact arise intrinsically from surface 
tension, being therefore unavoidable even in nominally 
"perfect" free-standing ferroelectric nanostructures. 



Model 

The development of the model focuses on the geome- 
try of the single crystal lamellae first reported in ref.0, 
an example of which is shown in Figure 1(a). Figure 
1 (b) shows a schematic cross-section of our model rep- 
resentation for such BaTiOs lamella, with the surface 
encapsulation layers (SL's) that are potentially respon- 
sible for the formation of 90°domains present and their 
obeyance to Kittel behaviour. The model comprises a 
tetragonal ferroelectric/ferroelastic layer (FL) of thick- 
ness Afl, exemplified by BaTiOs, sandwiched between 
two cubic paraelectric surface layers (SL) of thickness 
Asl- The surface layers provide the stress for the cre- 
ation of the ferroelectric- ferroelastic domains in the FL. 
The assumption of cubic surfaces is supported by the ob- 



servation that barium titanate nanoparticles have a core- 
shell structure with tetragonal interior and cubic surfaces 
lfj 17]. However, the cubic symmetry is not a key fea- 
ture of the model: for the SL to impose stress on the inte- 
rior, it is just sufficient that it does not undergo the phase 
transition at Tq- A non- ferroelectric capsule/matrix was 
also assumed in recent phase-field simulations of domain 
patterns in ferroelectric nanostructures [18j . 

In the present model, the SL thickness is an adjustable 
parameter to be determined from experimental measure- 
ments. The cartesian z-axis is oriented perpendicular to 
the plane of the lamella, and the x and y axes coincide 
with the crystal axes of BaTiOs as shown in Fig. QJi. The 
equilibrium lattice parameters of the FL are assumed to 
be those of the bulk ferroelectric tetragonal BaTiOs crys- 
tal and, in order to minimize depolarization fields, the 
polarization will tend to lie within the XY plane, point- 
ing parallel to either the X or Y directions. In terms of 
notation, [co,ao,ao] if P is parallel to x or [ao,co,ao] if 
P is parallel to y. The formation of the SL with different 
equilibrium lattice constants [bo, bo, bo] results in elastic 
stress, provided by lattice matching at the FL and SL 
interface. 

Consider to start with the simplest possibility when 
the matching stress is uniform, i.e., there are no stress 
gradients. Because SL is much thinner than FL, the lat- 
ter remains undeformed, keeping the equilibrium BaTiOs 
lattice parameters ([co,ao,ao] or [&o,co,&o])- Deforma- 
tion concerns only the SL which, because of the matching 
conditions, should conserve the same XY-plane parame- 
ters as the FL but can relax in the Z direction (s' zz = 0). 
The corresponding deformation energy is caused by the 
in-plane misfit strains of SL and contains the tension and 
shear parts [H, H EH, Q : 
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where the diagonal components of the plain strain tensor 
of SL (s f a ,s f c ,0) = (s f xx ,Sy y ,s f zz ) are expressed via the 
mismatch parameters: 
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G is the shear modulus, v is the Poisson ratio and the 
factor 2 corresponds to two (top and bottom) SLs on both 
sides of FL. We assume that G « 40 - 55 • 10 9 7Vm- 2 E3 



and v « 0.28 [24[ are approximately the same for SL and 
FL. 

Such system can be unstable towards formation 
of the experimentally observed periodic structure of 
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90°ferroelastic domains X/Y/X/Y shown in Fig. [TJ Usu- 
ally, such domains have 45°domain walls and "head-to- 
tail" polarization contacts to avoid the formation of de- 
polarization charge. The reason of the instability is that 
domain formation can reduce the elastic energy (pQ) of the 
SL by allowing the misfit stress to gradually relax inside 
the FL, within an interfacial region whose thickness is 
of the order of the domain width d. The domains will 
appear when the elastic energy imposed onto the SL by 
the ferroelectric is bigger than the energy required for the 
formation of domain walls. We emphasize that, although 
the net effect of ferroelastic twinning is to minimize the 
macroscopic deformation of the SL, this twinning neces- 
sarily requires a strain gradient near the interface. This 
is because, while afar away from the interface the lattice 
parameters of the FL are those of the bulk ferroelectric, 
at the interface they must become closer to cubic in or- 
der to match the SL. Accordingly, an inhomogeneously 
strained region must appear for which flexoelectric effects 
may be important ([H 0, H3, E^). The flexoelectric 



contribution, however, has been left out of this model for 
the sake of simplicity. 

The situation described here is in fact quite similar to 
the appearance of periodic 90°domains in epitaxial ferro- 
electric films strained by thick undeformable substrates 
Til . 2oL 22| in which the film/substrate mismatch strain 
relaxes in the near-surface layer of the film. The only dif- 
ference is that the contact elastic SL in our case is thin 
and its deformation should be considered self-consistently 
with that of FL. In a way, the model discussed here can 
be seen as a generalization of the substrate/film mod- 
els for the case of a substrate that has finite thickness 
compared with the film. 

The mechanism of nonuniform strain relaxation in the 
domain-populated FL gives rise to two new energy con- 
tributions: the near-surface deformation energy of the 
ferroelectric layer, Wfl, an d the energy of the domain 
walls Wdw- The energy balance between Wsl, Wfl and 
Wdw optimizes the domain period 2d and the matching 
plane lattice constants changing periodically at SL/FL 
interface from (c, a) at X-domain to (a, c) at Y-domain. 

Before proceeding to the derivation of the total energy 
of the system 



W = 2Wfl + 2Wsl + W dw 



(3) 



(the factors 2 correspond to two-side SL and two-side 
relaxation near-surface layers in FL), we assume that the 
optimal domain width d is thinner than FL but thicker 
than SL : 



We express the deformation energy Wfl of the ferro- 
electric layer in terms of periodically changing strains of 
domains (s a , s c )/(s c , s a ) at the surface of the FL, 
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These are taken as variational parameters in the gen- 
eral expression for the elastic energy of the 2-dimensional 
crystal, periodically strained as (s a , s c )/(s c , s a ) and with 
relative domain population equal to <j) (See also Eq.(30) 
in [H): 
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Adapting Eq.(|6]) to the case of equally-populated (cj) = 
1/2) and thin (d < Afl) domains, the universal dimen- 
sionless function can be simplified as: f a « ApL/^d^ 
7C(3)/8tt 3 Q, and therefore: 
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Physically, the first term in (0) corresponds to the av- 
erage interface-induced strain that propagates through 
the whole thickness Afl of the FL. The second term 
is produced by the superposition of alternative strains 
of an infinite series of domains, compensating inside the 
FL and relaxing in the near-surface gradient layer of a 
thickness of order d. Technically this sum is expressed 
via the zeta- function C(3), like the electrostatic energy of 
alternative depolarization charge compensation in Kittel 
formula for ferroelectric domains [H, [ij . 

Before minimizing ([7J) we can ensure first the vanishing 
of the largest first (volume) term, selecting s a + s c = 0; 
this is justified because the domain pattern can compen- 
sate for shear strains (s a — s c ) but not for volume changes. 
The number of variational parameters then reduces to 
one: s = s a = — s c and the second (surface) term takes 
the form: 



W FL = 2x~ 1 dGs 2 . 



(8) 



Consider now the deformation energy Wsl of the SL 
subjected to 2d— periodic domain-induced surface strains 
( s a5 s c)/( s cJ s o)j taking into account that s' a and s' c are 
expressed via variational parameters s a and s c §5§ as: 



A SL <d< A FL . 



(4) 



Observed in experiment and discussed in detail later, 
such hierarchy simplifies the calculation of the different 
contributions to ([3]). 
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If SL is thinner than the domain structure period: 
(Asl <^ d) the periodic surface strain does not man- 
age to relax across the SL. Then, SL can be presented 



as a piecewise (s' a , s' c ) / (s' c , s' a ) strained film, having uni- 
form deformation for each section. Summing the given 
by the first part of Eq. (pp) elastic energies from all the 
pieces (that are equal because of the s' a <-> s' c symmetry) 
and taking into account the discussed above constrain 
s a = — s c = s we present the elastic energy of SL as a 
superposition of tension and shear contributions: 



W SL = ^As L G^-^(5 ab 



S cb ) 2 + -A SL G(2s-5 ac ) 2 . 
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Note that the tension contribution: 
W T 
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does not depend on strain variational parameter s and 
coincides with the tension energy of the uniformly de- 
formed SL in (pQ). 

Combining (|5J) and (fTD|) with the energy of domain 
walls 
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(a is the surface energy density of the domain walls), the 
total energy (j3]) is: 



W = W T + Cr- 
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where r = 2s/S ca and 

d = ^kDA fl = 2.7y r DA 

FL- 

with the length-scale parameter 

D = ——5- « 0.5nm, 



(13) 
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(15) 



that can be interpreted as " domain wall half-thickness" . 
This value is ca. one-two lattice constants and agrees 
with previously reported theoretical and experimental 

values M Ei El IB. 

Minimization of (fT3|) over r and d gives the expression 
for optimal domain width 



d = 



do 



(16) 



When do <C mAsl, one can neglect the responsible 
for the upward curvature denominator in ([T6]) and ob- 
tain the Kittel-type dependence ([T4|) for cI(Afl)- This 
explains why our experimental results could be reason- 
ably well fitted assuming a simple square root depen- 
dence of domain periodicity on film thickness [7]. Note 
however that we are still in the thin surface layer limit 
Asl < do. The opposite (irrelevant for our system) limit 
do < Asl will lead to the Roytburd-Pompe-Pertsev situ- 
ation of thick substrate with domain pattern also obeying 
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FIG. 2: Different stages of a FIB process to fabricate a 
BaTiOs lamella for cross-sectional view: (a) schematic draw- 
ing of a lamella milled by FIB; prior to milling, a protective 
rim of either Au or Pt (pale blue strip in the top edge of the 
lamella in the sketch) was deposited in order to help mini- 
mize Ga damage during milling, (b) FIB image of a lamella 
tilted at 45 for better viewing, (c) drawing of the process to 
cross-section a lamella, and (d) cross-sectional TEM image of 
a lamella milled by FIB. 



the Kittel law but with a different numerical constant. 
When do — > kAsl^ the domain width diverges, implying 
a transition to a mono-domain state. Domains therefore 
exist only when: 



do < xA S L 
or, taking into account (jHj) and (fTBj) when: 

A 2 
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where D is estimated as (fT5|) . This represents a rather 
narrow constraint, providing assumption (|4]) is satisfied. 

Substitution of the optimal parameters r and d into 
(fT3|) gives the energy of the domain state 



W = W T + 2 
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which is smaller than ([T|) (since 2x(l — x/2) < 1 with 
x = do I h Asl) m the limits of applicability of the the- 
ory. This confirms the instability of the ferroelectric free- 
standing lamella with surface tension towards domain 
formation. 



Comparison with experiment 

An attempt was made to use the model developed 
above to describe the variation in domain periodicity ob- 
served for single crystal BaTi03 lamellae quantitatively. 
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FIG. 3: (a) HRTEM image of a cross-sectional BaTi0 3 
lamella after FIB milling showing a 20nm thick damaged layer 
at the Pt/BaTi03 interface (the Pt epilayer was deposited to 
preserve the original surface structure associated with FIB 
processing of the lamellae) (b) Zoom in on the damage layer, 
showing its amorphous structure. 



To do this it was noted that the previously published do- 
main period data [7] had all been taken from lamellae for 
which there had been no attempt to repair surface dam- 
age caused by focused ion beam (FIB) processing. It was 
expected that physical 'encapsulation layers' of amor- 
phous BaTi03 should exist on the top and bottom lamel- 
lar surfaces. To establish the thickness of the physically 
damaged layers, cross-sectional transmission electron mi- 
croscopy was used, with the lamellar cross-sections pre- 
pared by FIB according to the schematic shown in figure 
[21 As can be seen in figure [3l a surface layer of amorphous 
material, approximately 20nm in thickness, was indeed 
observed. In conjunction with energy dispersive x-ray 
data, this layer was categorized as a gallium-impregnated 
barium titanate glass. 

Substitution of this glassy layer thickness as that of 
an encapsulation layer (Asl) an d u sing a domain wall 
energy density of a =3xlO _3 Jm _1 [29[ produced a re- 
markably good quantitative description of the observed 
domain periodicity data, as can be seen in figure [H This 
strong agreement was obtained without any free fitting 
parameters, and one might naturally conclude that all of 
the encapsulation suffered by the BaTiC>3 lamellae was 
indeed due to the constraint from the ion beam damaged 
layers. 

However, we have spent some considerable effort to 
develop processing methodologies to repair the surface 
damage caused by focused ion beam processing. Ther- 
mal annealing in air at 700 °C has been seen to both re- 
crystallise the damage and expel the implanted gallium 
(forming thin gallium oxide platelets), recovering pris- 
tine single crystal BaTiOs [23, [HI (see figure [5]). If 
the thermal annealing is performed in oxygen, then func- 
tional measurements even suggest that the permittivity 
of the surface region is the same as that seen in bulk 
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The domain structure seen in a lamella which had been 
thermally annealed and is expected to have no surface 
damage is shown in figure [5]-c. The domain appearance 
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FIG. 4: (Colour online) Solid data in black colour: experi- 
mentally measured domain width as a function of thickness 
for free-standing single crystal films of ferroelectric BaTiOs 
with a dead layer thickness of ~ 20 nm. The red colour point 
represents the domain width after annealing at 700 °Cfor 1 h 
in air. For a lamella thickness of 150 nm, the domain width 
changes from 76.6 ±6.5 nm (before annealing) to 88.4 ±3.5 
nm (after annealing). Dotted line, calculation using eq. (fT6|) 
assuming Asl — 20nm and a domain wall energy density of 
a ~ 3 • 10 -3 Jm" 1 [29]. The agreement is good without free 
fitting parameters 



is almost identical to the unannealed sample in Figure 1, 
with the domain walls still of 110-type, indicating that 
90°domain sets have again formed. While gallium oxide 
platelets have precipitated on the annealed surface (the 
white blotches in figure [5]-c), these platelets certainly do 
not form a continuous layer, and it seems unlikely that 
they can provide the homogeneous stress needed to in- 
duce the 90°domain sets observed. Further, when the 
periodicity of the domains is analysed and compared to 
that obtained in unannealed lamellae, only a slight in- 
crease in domain period is observed; the increase in do- 
main period is in fact consistent with an effective increase 
in the thickness of the lamella by 40nm (that associated 
with the recrystallisation of the amorphous glassy regions 
described above). The essential physics at play therefore 
appears to be unchanged even when the extrinsic damage 
layer is removed. 

Overall, then, it appears that the quantitative agree- 
ment between the model and our data may be fortuitous, 
and that the stress responsible for the domains is not nec- 
essarily coming from the glassy barium titanate. Instead, 
a more fundamental source of stress must be at play for 
all of the FIBed single crystal barium titanate lamellae 
investigated to date. This suggests the existence of an 
intrinsic surface relaxation layer, probably due to surface 
tension. 

The existence of surface lay ers in BaTiOs has been 
known for quite some time [ll [13, 0, EE S, [H, 0, EI 
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with tetragonal SL provided that the tetragonal axis is 
out-of-plane. 

On the other hand, the elastic energy stored by the sur- 
face layer is proportional to its thickness. The outstand- 
ing question, then, is whether an intrinsic and very thin 
SL due to surface tension can lead to the same domain 
size as would be expected from the thicker encapsulation 
layers seen in our unnanealed samples. 

Indeed, EqUU states that the domain size is essentially 
do, independent of the SL thickness; however, this is 
equation is only valid when do <C kAsl-, which is not 
true if Asl ~ lnm. On the other hand, our model has 
implicitly assumed that the stiffness (the shear modulus 
G) is the same for FL and SL, and this is unlikely when 
the SL is an intrinsic surface-tension layer; for these, the 
bonds are known to be shorter [42[ and the SL should 
be expected to be harder than the FL. It is relatively 
straight-forward to incorporate the different shear mod- 
ulus of the surface layer (Gsl) and the ferroelectric layer 
(Gfl) onto the model, by simply substituting G for Gfl 
and Gsl in eqs. dHJ and (fTOj) respectively. Minimization 
of the total energy then leads to the new generalized ex- 
pression for domain size: 



d 



(c) 

FIG. 5: (a) HRTEM image of a cross-sectional BaTi0 3 
lamella after FIB milling and annealing at 700C in air. (b) An 
inverse FFT image of the Au/BaTi03 interface showing the 
complete reconstructed structure by post-annealing. Lattice 
fringes are evident up to the boundary with an Au protec- 
tive layer, with no evidence of amorphization (c) Plan-view 
TEM image of the annealed lamella, showing the expelled Ga 
as clusters on the surface and the persistence of the regular 
domain pattern. 



42l . |43[, yet there is surprisingly little agreement about 
their properties. Experimentally, their thickness seems 
to be on the region of 5-10 nm or more 16|, 17, 3(| 37, 38[ 
, whereas first principles calculations give a much smaller 
value, about lnm 

m Ed, IHIH • in some works the SL 

is found to be tetragonal at all temperatures, even above 
Tc [36[ [37j, where in others it is cubic even below Tc 
[Til, nq. In fact, the structure of the SL can be rather 
complex and depends on processing conditions . Nev- 
ertheless, what is important from the point of view of our 
model is not so much the exact symmetry of the surface 
, but the fact that it does not undergo the same ferro- 
electric/ferroelastic phase transition as the inside of the 
film. Thus, when the film becomes ferroelectric, it auto- 
matically becomes stressed by the untransformed surface 
layer. The present model requires only the induced stress 
is isotropic or orthotropic, and this is the case not only 
with cubic SL's, but also with either amorphous SL or 
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Gsl 



(20) 



This expression is almost identical to (fT5|) except for 
the appearance of the factor Gfl /Gsl in the denomi- 
nator. This factor can compensate for a reduced thick- 
ness of the surface layer insofar as its hardness is greater 
than that of the ferroelectric layer. Presently we have 
no quantitative estimates for the hardness of the intrin- 
sic epilayer, and we very much encourage the theoretical 
community to perform first principles calculations of the 
value of Gsl- What we can say is that, if the thickness of 
the SL is only l-2nm, as suggested by the ab-initio calcu- 
lations [39|, [13, 41, 42|, its shear modulus would need to 



be roughly 10 times bigger than that of the FL to have 
the same effect on domain size as our 20nm extrinsic en- 
capsulation layer. If, however, the true thickness of the 
intrinsic SL is ca. lOnm, as s ugg ested by experimental 
measurements 1(3, 17, 36, 37, 38j then the SL need not 
be more any more rigid than the FL. 



Conclusions 

In sum, we have shown that, even in the absence of 
rigid substrates or any other source of external stress, 
ferroelastic twinning can appear due to the self-stress im- 
posed by surface layers. The importance of such layers 
obviously increases as the size of the system decreases 
such that this effect becomes particularly important at 
the nanoscale. Furhermore, there need not be an extrin- 
sic surface layer for the ferroelastic domains to appear; 
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surface tension, which is intrinsic and therefore unavoid- 
able, can also provide the necessary stress for domain for- 
mation. Finally, we note that epilayers can be expected 
to be important not only for isolated nano-ferroelectrics, 
but also in macroscopic devices such as ceramic capaci- 
tors made with nano-powders or core-shell grains. 

This work was supported by the EC project FP6- 
STREP-MULTICERAL and by the French-UK collab- 
oration program "Alliance" . 
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